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Abstract 

A strictly truncated (weak-coupling) perturbation theory is applied to the 

attractive Holstein and Hubbard models in infinite dimensions. These results 

are qualified by comparison with essentially exact Monte Carlo results. The 

second order iterated perturbation theory is shown to be quite accurate in 

calculating transition temperatures for retarded interactions, but is not as 

accurate for the self energy or the irreducible vertex functions themselves. 

Iterated perturbation theory is carried out thru fourth order for the Hubbard 

model. The self energy is quite accurately reproduced by the theory, but the 

vertex functions are not. Anomalous behavior occurs near half filling because 

the iterated perturbation theory is not a conserving approximation. 
Pacs:74.20.-z, 71.27.+a, and 71.38.+i 
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I. INTRODUCTION 

Recently Metzner and VollhardtS showed that the fermionic many-body problem simpli- 
fies in the limit of infinite spatial dimensions. The d ^ oo limit is taken in such a fashion 
that the dynamics are local in space, and the lattice problem can be mapped onto a self- 
consistently embedded Anderson impurity model. The self energy and irreducible vertex 
functions of the impurity problem are employed to determine the Green's functions and 
susceptibilities of the infinite-dimensional lattice. This many-body problem can essentially 
be solved exactly by using the quantum Monte Carlo (QMC) techniques of Hirsch and FyJ 
to extract the self energy and irreducible vertex functions of the relevant impurity problem. 
A variety of models have already been examined in this fashion: the Hubbard modelH^; 
periodic Anderson mo del0; and the Holstein modelu. 

The infinite-dimensional limit provides a unique testing ground for various approximation 
techniques, since one can compare them to the benchmark QMC results. Previous work has 
concentrated on the strong-coupling limiti§ (perturbation theory in the kinetic energy) 
which can be shown to be an excellent approximation in the region of moderate to strong 
coupling. 

Calculations in the weak-coupling limit are not under as good control. The main em- 
phasis so far has been on conserving approximationslliFii (such as the fiuctuation-exchange 
approximation) to the many-body problem. Conserving approximations are popular because 
they satisfy the requisite conservation laws for electronic charge, energy, and momentumS. 
It turns out, however, that conserving approximations suffer from some serious fiaws: (1) 
truncated approximations (thru fourth order) do not show a turnover in as the strong- 
coupling regime is approached; (2) infinite summations of certain classes of diagrams can 
produce a turnover in Tc, but the self energy is overestimated, and many-body features 
(such as the upper and lower Hubbard bands) do not appear in the single-particle spectral 
function. Similar deficiencies where found by White in his examination of the Anderson 
impurity model itsel£^. 
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An alternate approximation scheme, that does not suffer from the above problems (but 
is not conserving), has been introduced by Georges and Kotliar§ and is known as iterated 
perturbation theory (IPT). In this case, the perturbation theory is strictly truncated at 
a finite order (as opposed to a conserving approximation that sums the infinite class of 
diagrams corresponding to self energy insertions into the dressed Green's function). The 
motivation for this approximation comes from the pioneering work of Yosida and YamadaS, 
who found that truncated perturbation theory was an extremely accurate approximation for 
the Anderson impurity model because the sum total of all classes of diagrams at a given 
order is an order of magnitude smaller than any of its constituent parts. This remarkable 
near cancellation of diagrams was further verified to high order by Zlatic and Horvaticlli who 
found that the exact Bethe ansatz solutions can be rearranged into a power series whose 
coefficients die off rapidly. 

Horvatic and ZlaticS generalized Yosida and Yamada's analysis off of half filling and 
found that the strictly truncated perturbation series was also accurate in the doped case. 
However, anomalous behavior begins to enter if the coupling strength is increased too far 
and the electron concentration is close to half filling. Ferrer, Martin-Rodero, and Flore£ll 
realized that the reason why the perturbation theory fails more rapidly away from half 
filling is because the second-order perturbation theory correctly reproduces the atomic limit 
at half-filling, but not away from half filling. They proposed an ad hoc interpolation scheme 
that properly reproduces the atomic limit everywhere and removes some of the anomalous 
behavior of the Horvatic and Zlatic approximation. An alternate interpolation scheme has 
been proposed by Nealil based upon functional integral techniques which appears to be 
quite accurate for all values of the interaction strength but is restricted to half filling. 

In this contribution, we will examine the iterated perturbation theory through second 
order for the Holstein model and through fourth order for the negative-U Hubbard model. 
Comparison of the self energy, the irreducible vertex functions, and the transition temper- 
atures will be made to the exact QMC solutions. In Section II the infinite-dimensional 
formalism is reviewed along with the details of the IPT approximation. Section HI contains 



the results for the Holstein model, while Section IV holds the results for the Hubbard model. 
Conclusions and discussion follow in Section V. 



II. FORMALISM 

The electron-phonon interaction is generally believed to be responsible for the supercon- 
ductivity in most low-temperature superconductorJ^lii. The conventional theory used to 
describe superconductivity was introduced by Migdalil and Eliashbergil. Migdal-Eliashberg 
(ME) theory ignores vertex corrections which can be shown to be proportional to powers of 
the Debye frequency divided by the Fermi energy, and hence are small corrections for most 
real materials. 

Here, iterated perturbation theory is examined in a systematic fashion for the attractive 
electron-phonon interaction described by the Holsteinll model (the effects of Coulomb re- 
pulsion are explicitly neglected). The Holstein model consists of conduction electrons that 
interact with localized (Einstein) phonons: 

^ = - E (SvCfco + cLc,.) + Y.(9Xj - /i)(n,T + - 1) + ^Mfi^ Y,:,^. + W?± 

(1) 

where c]^ {cja) creates (destroys) an electron at site j with spin a, rij^ = (^j„Cj^ is the electron 
number operator, and Xj (pj) is the phonon coordinate (momentum) at site j. The hopping 
matrix elements connect the nearest neighbors of a hypercubic lattice in (i-dimensions. The 
unit of energy is chosen to be the rescaled matrix element t*. The phonon has a mass M 
(chosen to be M = 1), a frequency Q, and a spring constant k, = MQ"^ associated with it. 
The electron-phonon coupling constant (deformation potential) is denoted by g; the effective 
electron-electron interaction strength is then the bipolaron binding energy 

The chemical potential is denoted by fi and particle-hole symmetry occurs for /i = 0. 



In the instantaneous limit where U remains finite and g and fl are large compared to the 
bandwidth {g,Q oo,U = finite), the Holstein model maps onto the attractive Hubbard 
modelS 

t* 11 

H = -—7^ E + cLcj.) - /iEKt + ^ii) + u 5IKt - - 7^) (3) 

with U defined by Eq. (g). 

The infinite-dimensional limit of Metzner and VollhardtS is taken {d —>■ 00), in which 
the electronic many-body problem becomes a local (impurity) problem that retains its com- 
plicated dynamics in time. The hopping integral is scaled to zero in such a fashion that 
the free-electron kinetic energy remains finite while the self energy for the single-particle 
Green's function and the irreducible vertex functions have no momentum dependence and 
are functionals of the local Green's function0'@@. This limit retains the strong-correlation 
effects that arise from trying to simultaneously minimize both the kinetic energy and the 
potential energy. 

The many-body problem is solved by mapping it onto an auxiliary impurity problemiE 
in a time-dependent field that mimics the hopping of an electron onto a site at time r and 
off the site at a time r'. The action for the impurity problem is found by integrating out 
the degrees of freedom associated with other lattice sites in a path-integral formalism. § The 
result is an effective action 

Seff. = /o^ dr dT'cUT)G,\T - t'Mt') + /cf dT[gx{r) - [4(r)c.(r) - i] 

+ IM drin^x^r) + x\t)] (4) 

where Gq^ is the "bare" Green's function that contains all of the dynamical information of 
the other sites of the lattice. The interacting Green's function, defined to be 

is determined by Dyson's equation 

G-' = G-\iUr,) = Go i(z^„) - S(2a;„)- (6) 



A self-consistency relation is required in order to determine the bare Green's function 
Gq. This is achieved by mapping the impurity problem onto the infinite-dimensional lattice 
thereby equating the full Green's function for the impurity problem with the local Green's 
function for the lattice 

Gjj{iUn) = G(k, iuri) = Ylli^^n + - Eik) - = F^[iuJn + - (7) 

k k 

Here Foo{z) is the scaled complimentary error function of a complex argument 

]^ roo expf y'^) 2 

Foo{z) = —j= / dy = — isgn[Im(2;)]-\/7re~^ erfc{— zsgn[Im(2)]z}. (8) 

a/tt j-oo z — y 

The dynamics of the (local) impurity problem are identical to the dynamics of the Anderson 
impurity modeliiii3ii0i. This many-body problem can either be solved exactly with the 
QMC algorithm of Hirsch and FyJ, or it can be approximately solved by employing a 
truncated perturbative expansion. The impurity is self-consistently embedded in the host, 
since it must satisfy the self-consistency relation in Eq. (^. 

It is important to note that since one does not a priori know the bare Green's function 
Go^ in Eq. (§, one must iterate to determine a self-consistent solution for the Green's 
function of the infinite-dimensional lattice. This is achieved by calculating the self energy 
as a functional of the bare Green's function Gq, and then determining the new local Green's 
function from the approximate self energy and Eq. (|^). The new bare Green's function is then 
calculated from the Dyson equation in Eq. @. This process is iterated until convergence is 
achieved [for the perturbation theory approximation the maximum variation of each G{iuJn) 
is less than one part in 10^ which typically takes between 5 and 30 iterations, for the QMC 
calculations the algorithm is iterated from 7 to 9 times]. Note that the iteration required 
to determine a self-consistent bare Green's function Go should not be confused with the 
the iterative techniques used to self-consistently sum all of the self-energy insertions in a 
self-consistent perturbation theory — the perturbation series is always strictly truncated at a 
finite order here. 

The QMC algorithm proceeds by discretizing the imaginary-time interval from to /3 
into L time slices of equal width At = f3/L and evaluating the relevant path integrals in 
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a grand canonical scheme. Both local moves, in which the phonon coordinate is shifted by 
a different amount at each time slice, and global moves, in which the phonon coordinate is 
shifted by a uniform amount for every time slice, are incorporated into the QMC algorithmii. 
The values of L used ranged from 20 to 160 with the largest values of Ar reserved for the 
lowest temperatures (usually Ar was fixed at 0.4). No sign problem was found at any filling. 

Static two-particle properties are also easily calculated since the irreducible vertex func- 
tion is local0. The static susceptibility for CDW order is given by 



CDW 

A. 



;q) = ^ E e^'^-(^^-^'^-)T fdr f dr'[{n,^{r)nUr')) - Mr))Mr'))] 



mn 



at each ordering wavevector q. Dyson's equation for the two-particle Green's function 
becomeJiS 

Xmn (q) ~ Xm(q)^mri. ~ T^^^ Xm(q)-^mp Xpn (q) ' (lO) 

P 

with r^^'^ the (local) irreducible vertex function in the CDW channel. 

The bare CDW susceptibility Xn(q) iii Eq. (|10|) is defined in terms of the dressed single- 
particle Green's function 



Xn(q) = -^EGn(k)G„(k + q) 



1 1 , e-y" 



dy- 



V^i/l -X2(q) J-^ iuJn + fi - - y 



iujri + j U- S„ -X(q)y 
1-X2(q) 



and all of the wavevector dependence is included in the scalaicilO X(q) = cosq^/d. 
The mapping q ^ X{(\) is a many-to-one mapping that determines an equivalence class of 
wavevectors in the Brillouin zone. "General" wavevectors are all mapped to X = since 
cosqj can be thought of as a random number between —1 and 1 for "general" points in 
the Brillouin zone. Furthermore, all possible values of X (—1 < X < 1) can be labeled by 
a wavevector that lies on the diagonal of the first Brillouin zone extending from the zone 
center {X = 1) to the zone corner (X = — 1). The presence of incommensurate order in the 



attractive Holstein model is restricted to a very narrow region of parameter spaceuiB so only 
the "antiferromagnetic" point X = — 1 is considered for CDW order. The integral for x^(X) 
in Eq. (|Il|) can then be performed analytically^ Xni-^ = = —Gn/ii^n + — The 
irreducible vertex function T^^^ is either directly calculated in a perturbative expansion 
(IPT) or is determined by inverting the Dyson equation in Eq. (|I^) (QMC). 

A similar procedure is used for the singlet s-wave SC channel. The corresponding defi- 
nitions are as follows: The static susceptibility in the superconducting channel is defined to 
be 

= ^ E e^'^-^^^^^'^^^rrfr rrfr'(c,Kr)c,i(r)4(r')cI,(r')) 

iV-j-jT-j J J 

ti-j—r^k 

^Tj2x'^{q,tu;r,,,tu;r^)=TY.X'l'ni^) , (12) 

mn mn 

for superconducting pairs that carry momentum q; Dyson's equation becomes 

X^n(q) = xl'icD6^n - T5:x^'(q)r^pXpf (q) , (13) 

p 

with r,;^*^ the corresponding irreducible vertex function for the SC channel; the bare pair-field 
susceptibility becomes 



Xn(q) = ^EG'n(k)G_„_i(-k + q) 



1 



/oo 
dy 
-oo 



^ Jl - X2(q) iuJn + fJ' - - y 



(14) 

X2(q 

with the special value Xn'(X = 1) = — ImG'„/Im(ici;„ — S„) for the SC pair that carries no net 
momentum; and finally the irreducible vertex function is also either directly calculated in a 
perturbative expansion (IPT) or is determined by inverting the Dyson equation in Eq. (p!3D 
(QMC). 

There are two different approximation schemes that can be used for determining the 
irreducible vertex functions in a perturbative expansion. The first scheme, denoted IPT, 
is in the spirit of Yosida and Yamada's original analysiJi^: the vertex functions are also 
expanded in a perturbation series that is strictly truncated at a finite order. In the second 
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scheme, denoted IPT*, the fully dressed local Green's functions are used to calculate the 
irreducible vertex functions. In general, one might expect the IPT approximation to be 
more accurate than the IPT* approximation from the work of Yosida and YamadaS, but 
the irreducible vertex functions were not investigated in their analysis of the Anderson 
impurity model, so it is not a priori clear which will be a better approximation. 

The transition temperature of the infinite-dimensional Holstein model is now found by 
calculating the temperature at which the relevant susceptibility diverges (CDW or SC). One 
can determine this transition temperature by finding the temperature where the scattering 
matrix (in the relevant channel) 

has unit eigenvalueS. Note that the fully dressed Green's functions are always used in 
calculating the bare susceptibility x° because the bare susceptibility is a property of the 
infinite-dimensional lattice, and not the Anderson impurity problem, and we expect to find 
better agreement with the QMC results if we use the correct bare susceptibility. 



III. HOLSTEIN MODEL 

The diagrammatic expansion for the IPT is depicted in Figures 1 and 2. The wiggly 
lines denote the bare phonon propagator, and the straight lines denote the bare Green's 
function Go- Figure 1 shows the self energy for the IPT approximation through second 
order. The self energy includes, respectively, the Hartree term (which is a constant that is 
reabsorbed into the chemical potential), the Fock term, the second-order term that dresses 
the phonon propagator, the lowest-order vertex correction, and the second-order rainbow 
diagram (which corresponds to the self-energy insertion of the first-order Fock diagram into 
the first-order Fock diagram). 

To be more explicit, the self energy for the second-order IPT approximation is 

Sn = -^T^ ^' Go(^^.) 
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n—r 

^u^rY^^^^^^fAM'-.) (16) 

which includes the Fock diagram contribution and the three second-order contributions of 
Figure 1. The bosonic Matsubara frequency vi is defined to be vi = 21ttT. 

The self-consistency step involves first determining a new local Green's function G„ from 
the integral relation in Eq. (|^) with the approximate self-energy of Eq. (TB). Next the 



new bare Green's function is found from the Dyson equation in Eq. (|^) using the same 
approximate self energy. A new self energy is then calculated from the new bare Green's 
function [using Eq. ([T6|) ], and the process is iterated. This iteration process terminates when 
the maximum deviation in the local Green's function is less than one part in 10^. 

Once the Green's functions and self energies have been determined, the irreducible vertex 
functions can be calculated for the CDW or SC channels. In the CDW channel [see Figure 
2 (a)] one must include both direct and exchange diagrams as well as the vertex corrections. 
The result is 



fi2 + zy2 



, ^ , , , , fi2 ^2 ^2 

^7^^7—2 2U^TJ2Go{^uJr)Goi^oo^-r 



^2 ^2 

for the IPT approximation. The IPT* approximation has the same functional form as in 
Eq. dT^), but the bare Green's function Go is replaced by the dressed Green's function 
G (this functional form is the same as the one used for the conserving approximation0) . 
Note that the vertex corrections (arising from the first-order exchange diagrams) modify the 
interaction in the CDW channel so that it properly interpolates between the zero frequency 
limit r*^"^^ — > 2f/ and the infinite frequency limit r*^^^ U. At an intermediate frequency, 
the CDW interaction strength has a complicated temperature dependence. In the SC channel 
[see Figure 2 (b)] one finds 
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Q2 

r^^ = u - 

it -r i^rn-n 

Vt^ Vt^ Vt^ 

+ ^7^T^Go(^^.)go(^^n^-n+.) 2 , ^2 [^0^-72 02 , ^2 " Q?-T2 

- f/^T^Go(^^^r)G'o(2U;m+n+r+l) 



m+r+1 ^' "T '^n+r+l 



for the IPT approximation. Once again, the IPT* approximation has the same functional 



form as in Eq. (18), but the bare Green's function Go is replaced by the dressed Green's 
function G. 

At half filling the Holstein model interaction is particle-hole symmetric, so the Green's 
function and self energy are purely imaginary and the vertices are real. The self energy 
can be expressed by Ti{iuj„) = iuJnZ{iuJn), with Z{iuJn) the renormalization function for 
the self energy. At half filling, both Xmi-^ — ^"^^ ^^e maximal eigenvector of the 
scattering matrix are also even functions of Matsubara frequency, so the only contribution 
of the irreducible vertex function to the eigenvalue of the scattering matrix comes from the 
even Matsubara frequency component [rm,n + r-m-i,n]/2. The same result holds for the SC 
channel at all fillings, because both the maximal eigenvector and the bare susceptibility are 
symmetric functions of Matsubara frequency. 

A comparison of these approximations to the exact QMC results for the electronic self 
energy and one column of the irreducible vertex function in the CDW and SC channels is 
made in Figures 3-5 for three different interaction strengths at half filling. The phonon 
frequency is set to be approximately one-eighth of the effective bandwidth {Q/t* = 0.5) as 
was done in the QMC solutions!. The energy cutoff is set to include 256 positive Matsubara 
frequencies for the perturbative approximations. 

At weak coupling {g = 0.4, T = 0.0625, Figure 3), the second-order IPT approximation 
is a reasonable approximation to the electron self energy [Figure 3 (a)] and to the CDW 
irreducible vertex function [Figure 3(b)]. Note that the IPT underestimates the self energy 
(an effect that raises Tc) and overestimates the magnitude of the vertex (an effect that also 
raises Tc), so one expects it to overestimate the CDW transition temperature (both the 
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IPT and the IPT* approximations yield similar results for the irreducible vertex function). 
Clearly third-order diagrams play an important role even at this weak a value of the coupling 
strength. In this sense, the IPT is less accurate for retarded interactions than it is for 
instantaneous interactions (see Section IV). 

As the coupling strength is increased to the point where a double-well structure begins 
to develop in the effective phonon potential of the QMC simulations of Ref. i| ((? = 0.5, 
T = 0.111) one can see strong-coupling effects begin to enter into the QMC vertex. The 
IPT still underestimates the self energy [Figure 4 (a)], but as the strong-coupling regime 
is approached, the vertices become increasingly attractive [Figure 4 (b)]. This enhanced 
attraction is not represented by either the IPT or the IPT* approximations. Note that in this 
regime the self energy and the vertices are underestimated, which implies that the transition 
temperature may be determined more accurately than expected, because these effects tend 
to cancel each other out in determining Tc. The underestimation of the magnitude of the 
vertex will eventually cause to drop due to the large self-energy renormalization. 

When one is well into the strong-coupling regime {g = 0.625, T = 0.167, Figure 5) the 
deficiencies of the IPT approximation become more apparent. The self energy is underes- 
timated by almost a factor of two [Figure 5 (a)], and the vertices are underestimated by 
almost a factor of three at intermediate frequencies [Figure 5 (b)]. Even at this large a value 
of the coupling, both the IPT and IPT* approximations remain similar to each other. At 
this point one would say the approximate theory is failing to faithfully represent the exact 
solution. 

A comparison of the IPT approximation to either a second-order conserving approxima- 
tion or to ME theory0 shows that the IPT is a superior approximation for both the self 
energy and the vertices, but none of these approximate methods is accurate over a wide range 
of interaction strength, indicating once again the importance of the third-order diagrams. 

At half filling, the Holstein model always has a transition to a CDW-ordered phase at 
q = (tt, tt, tt...) (X = — 1). The transition temperature to this commensurate CDW is plotted 
in Figure 6 as a function of the interaction strength. Both second-order IPT approximations 
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are compared to the QMC simulationsi. The IPT and the IPT* approximations are very 
accurate in determining T^; the peak position and height are both reproduced well. This 
result is definitely a numerical coincidence, since the self energy and vertices are poorly 
approximated at the peak of the curve {g = 0.625). Furthermore, it is the self energy 
renormalization that ultimately causes to turnover, but drops far too rapidly in this 
strong-coupling regime. This result is similar to what was found for a strong-coupling 
approximation!: there the transition temperature dropped too rapidly in the weak-coupling 
regime. 

The IPT and IPT* approximations are superior to either a second-order conserving 
approximation or to ME theory in both a qualitative and a quantitative determination of 
the CDW transition temperature at half filling. They also have the correct limiting behavior 
as [/ — because they are second-order perturbation theories000. 

As the system is doped away from half-filling, the CDW instability remains locked at the 
commensurate point {X = —1) until it gives way to a SC instability (incommensurate order 
may appear in a very narrow region of phase space near the CDW-SC phase boundaryili 
but is neglected here). In Figure 7, the phase diagram of the Holstein model is plotted for 
three values of the interaction strength {g = 0.4, g = 0.5, and g = 0.625). The weak-coupling 
QMC datai (g = 0.4) are reproduced most accurately by the perturbative approximations, 
as expected from the comparison of the self energy and the vertices in Figure 3. As the 
coupling increases to g = 0.5 (where the double- well structure begins to form in the effective 
phonon potential!), the approximations become less accurate. The anomalous behavior of 
the CDW transition temperature not being a maximum at half filling already appears for 
the IPT approximation, but the IPT approximation is more accurate at determining the 
CDW-SC phase boundary. This anomalous behavior in the CDW transition temperature 
most likely occurs because the approximation is not a conserving approximation. Note that 
the SC transition temperatures are still accurately reproduced in this regime. 

When one is well into the strong-coupling regime {g = 0.625) both the IPT and IPT* ap- 
proximations have large anomalous behavior near half filling. The CDW-SC phase boundary 
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is also poorly approximated. Clearly both approximations are failing at this large a value 
of the coupling. 

Both the IPT and IPT* approximations yield virtually identical results for the Holstein 
model at a moderate phonon frequency. The self energy and the irreducible vertex functions 
are poorly reproduced, implying that third-order diagrams are important when the inter- 
action is retarded, but remarkably, the transition temperatures are reproduced quite well. 
At half filling both approximations show a peak in the CDW transition temperature. Off 
of half filling, the SC transition temperature is determined more accurately than the CDW 
transition temperature (indeed anomalous behavior enters into the CDW channel near half 
filling when the interaction strength is large enough). The phase boundary between CDW 
and SC order is also accurately reproduced. 



IV. HUBBARD MODEL 

The Hubbard model in Eq. (|^) is the infinite-frequency limit [Q oo) of the Holstein 
model. The Hubbard model has an electron-electron interaction that only occurs between 
electrons with opposite spins. This happens because of the cancellation of the direct and 
exchange diagrams which causes all electron-electron interactions between like-spin particles 
to vanish. The perturbation theory becomes much simpler in the Hubbard model case, 
because of this reduction of diagrams, and can be performed to higher order. Here the 
truncated IPT approximation will be carried out to fourth order, and will be compared to the 
QMC simulations! to determine their accuracy. Previous work has concentrated on second- 
order conserving approximation^lllli, third- and fourth-order conserving approximations^, 
or the fiuctuation-exchange (FLEX) approximation00. 

One expects that a truncated approximation will be superior to an infinite summation of 
random-phase approximation bubbles and particle-hole and particle-particle ladders because 
the many-body problem reduces to a self-consistently embedded Anderson impurity model, 
and the analysis of Yosida and Yamadalil has shown that the total fourth-order corrections 
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to the self energy are an order of magnitude smaller and opposite in sign to the fourth-order 
contribution of the FLEX approximation. The irreducible vertex functions should have 
similar effects, but have not yet been analyzed in detail. 

The diagrammatic expansion for self energy (in the IPT approximation) of the Hubbard 
model is given in Figure 8. Here, the solid lines denote the bare Green's function Gq and 
the dotted lines denote the instantaneous Coulomb interaction U. The first line includes 
the first-order Hartree contribution (which is a constant that is absorbed into a renormal- 
ized chemical potential), the second-order contribution and the third-order particle- hole and 
particle-particle ladders. The second line contains the fourth-order contributions from the 
RPA bubbles and the particle-hole and particle-particle ladders. The third and fourth lines 
include all of the remaining diagrams that enter at fourth-order. The fifth line contains the 
insertion of the second-order self energy into the second-order diagrams. The irreducible ver- 
tex functions are too cumbersome to represent diagrammatically. Explicit formulas for the 
electronic self energy of the Hubbard model thru fourth order and for the CDW vertex have 
been given before^ (except for the second-order self-energy insertions into the second-order 
diagrams which are easily determined). Once again IPT denotes the approximation where 
the irreducible vertex is is strictly truncated at a given order, while the IPT* approximation 
determines the vertex functions by using the dressed Green's functions. 

Since the Hubbard-model interaction is particle-hole symmetric, the half-filled band cor- 
responds to = 0, and the Green's functions are purely imaginary. The odd-order contribu- 
tions to the self energy all vanish and each of the fourth-order contributions on a given line 
in Figure 8 are identicalEl. Since the self-energy is an even function of U at half-filling, but 
the irreducible vertex function contains both even and odd powers of U, the only difference 
between a truncated approximation of order 2n and of order 2n + 1 is that the irreducible 
vertex function is larger for the odd-order approximation. Therefore, we expect that an even 
order approximation will underestimate the transition temperature (in weak-coupling) and 
an odd-order approximation will overestimate T^. 

One of the most remarkable results of Yosida and Yamada is that the total fourth-order 
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contribution to the self energy is an order of magnitude smaller than any of its component 
pieces. One can ask if the same near cancellation holds in the mapping to the infinite- 
dimensional lattice. The answer is affirmative, as can be seen in Figure 9, where each of 
the individual contributions to the self energy are plotted at half-fiUing for [/ — —1 and 
T — 0.05. The dashed line shows the second-order contribution to the self energy. The 
dotted lines are each of the four different fourth-order contributions. FLEX denotes the first 
fourth-order contribution corresponding to the fourth-order bubbles and ladders. Line 2 and 
line 3 denote the contributions arising from the next two lines respectively in Figure 8, and 
IPT is used to denote the contributions arising from the second-order self-energy insertions 
into the second-order diagram (the last line of Figure 8) . The total of all of the fourth-order 
diagrams is plotted with the solid line. One can see that the remarkable near cancellation of 
the fourth-order diagrams still holds in infinite- dimensions, and in this case the fourth-order 
contribution is the same sign, but more than an order of magnitude smaller than the FLEX 
contribution at fourth order. 

A comparison of the different approximation schemes versus the exact QMC results are 
given in Figures 10-12 for three different values of U. The second-order and third-order 
approximations employ an energy cutoff of 256 positive Matsubara frequencies; the fourth- 
order approximation uses 64 positive Matsubara frequencies. In Figure 10 (a) the self-energy 
renormalization function is plotted for the two different approximations at C/ — —1 {T — 
0.05). Note that the fourth-order approximation virtually reproduces the QMC results. In 
Figure 10 (b) the even component of one row of the irreducible vertex function for the CDW 
channel are compared for U = — 1 in the IPT approximation. All of the approximations 
are in reasonable agreement with the QMC results, but the accuracy is significantly reduced 
relative to what was found for the self energy. The IPT* approximation is compared to 
the QMC results in Figure 10 (c). These results are similar to the IPT approximation, but 
perhaps not as accurate. 

As the coupling strength is increased to U = —2 (T = 0.125) the self energy is becom- 
ing bracketed by the two low-order approximations, with the fourth-order approximation 
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overestimating the self energy, and the second-order approximation underestimating the self 
energy, as shown in Figure 11 (a). The CDW vertex is plotted in Figure 11 (b) for the 
IPT approximation. One can see strong-coupling effects beginning to enter as the QMC 
vertex becomes more attractive at low frequency transfer than predicted by the even-order 
approximations. The third-order IPT approximation is producing the correct qualitative 
shape of the vertex, but it is overestimating the magnitude by a large amount. The IPT* 
approximation is shown in Figure 11 (c). Here the two different methods for calculating 
the vertex yield quite different results. In the IPT* approximation, all orders are showing a 
weakening of the vertex at lower frequencies, which is not present in the QMC data. 

Finally, at a rather strong value of the coupling {U = —3, T = 0.167), where the tran- 
sition temperature reaches its maximal value, the self energy is still reproduced quite accu- 
rately by both orders of the IPT [see Figure 12 (a)]. The vertices are, however, quite poorly 
approximated. In Figure 12 (b) and 12 (c) the results for the IPT and IPT* respectively 
are presented. The QMC vertex is becoming strongly attractive here. Only the third-order 
IPT approximation has the correct qualitative behavior, but it overestimates the magnitude 
of the vertex. The IPT* approximation produces a fairly fiat frequency dependence to the 
vertex which is quite accurate at all but the lowest frequency transfers (where it is almost 
an order of magnitude too small). 

The conclusion that we can draw from this comparison of the self energy and irreducible 
vertex functions is that the IPT approximation is quite accurate for the self energy, but 
is a rather poor approximation for the irreducible vertex function. It is known that the 
reason why the IPT is so accurate for the self energy arises from the fact that it properly 
reproduces the atomic limit0'@, but it is possible that it does not produce the atomic limit 
for the irreducible vertex functions. If this were the case, one could try another ad hoc 
interpolation scheme that interpolates between the perturbative result for the vertex, and 
the atomic limit, to provide a more accurate approximation. 

Note also that since the self energy is reproduced quite accurately by the IPT, but the 
magnitude of the vertex is grossly underestimated, the transition temperature will drop too 
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rapidly in the strong-coupling regime (with the exception of the third-order IPT approxima- 
tion where the transition temperature will grow too rapidly). This is illustrated by plotting 
the transition temperature for the CDW transition at half-filling in the attractive Hubbard 
model as a function of U (Figure 13) and comparing to the QMC solution!. In the IPT 
approximation [Fig. 13 (a)], the peak in as a function of U is properly reproduced for 
the even-order approximations, but not for the odd-order one. The transition temperature 
decreases too rapidly in the strong-coupling regime because of the underestimation of the 
vertex. In the IPT* approximation [Fig. 13 (b)], all orders have a peak in Tc as a function 
of U. The third-order approximation is the most accurate here, because it has the most 
attractive vertex, but even in this case the quantitative agreement with the QMC is not as 
good as was found for the truncated conserving approximation^, where the transition tem- 
perature did not have a peak, but the underestimation of the self-energy compensated for 
the underestimation of the vertex to produce reasonably accurate transition temperatures 
into the strong-coupling regime. 

The doping dependence of the SC transition temperature is plotted in Figure 14 for 
both the second and third order IPT and IPT* approximations and U = —1. Even at this 
weak a value of the coupling strength, the second-order IPT and IPT* approximations show 
the anomalous behavior of an increase in Tc for small dopings off of half filling. A third- 
order approximation works much better (because third-order contributions to the self energy 
enter when the electron concentration is not equal to 1), and does not display the anomalous 
behavior (at this value of U). Once again, the IPT* approximation is more robust against 
the anomaly in T^, but all IPT approximations quickly fail as the coupling strength increases 
beyond U = —1 (because they fail at half filling as shown in Fig. 13). 



V. CONCLUSIONS 

The iterated perturbation theory approximation^ has been applied to the attractive 
Holstein and Hubbard models in infinite dimensions. This weak coupling theory is not a 
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conserving approximation because it corresponds to a strictly truncated perturbation series 
for the self energy. The vertex functions are either approximated with a strictly truncated 
perturbation series too (IPT), or with a truncated perturbation series that includes all self 
energy insertions (IPT*). 

At half filling, these approximate theories are accurate in reproducing the electronic self 
energy (for the Holstein model, the accuracy increases as the phonon frequency increases) 
but do not reproduce the correct behavior of the irreducible vertex functions in the strong- 
coupling regime. 

In particular, for the half-filled Hubbard model, the IPT is a very accurate approximation 
for the self energy of the Hubbard model, but is a poor approximation for the irreducible 
vertex functions a low frequency transfer. As a result, the transition temperature does in 
general display a peak (because of the large self-energy renormalizations) , but decreases too 
rapidly in the strong-coupling regime. 

The IPT is a much less accurate approximation for the self energy of the Holstein model 
at finite phonon frequency (indicating the importance of higher-order diagrams). The IPT 
is also a poor approximation for the vertex functions, but, remarkably, the transition tem- 
peratures are quite accurately approximated both at half filling and off of half filling (for 
small enough coupling strength). The peak position and height of Tc at half filling are re- 
produced well, as is the doping dependence of Tc up to the region of parameter space where 
a double-well structure begins to develop in the effective phonon potential. 

As both of these systems are doped away from half filling, the transition temperatures (to 
either a CDW or a SC) display the anomalous behavior of initially increasing with doping. 
This anomalous behavior first enters at some critical value of the coupling, and becomes 
worse as the coupling strength is increased. For the Holstein model (with fl = 0.5) this 
critical coupling lies near g 0.5, whereas for the Hubbard model it occurs at \U\ ~ 1. 
Third-order diagrams tend to push this anomaly out to larger values of interaction strength 
as expected because the third-order contributions to the self energy can be large as the 
system is doped off of half filling. 
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As a final summary, the IPT is a highly accurate approximation for the self energy when 
the interaction is instantaneous, but generally underestimates the magnitude of the vertex 
functions at low frequency transfer; it does not approximate the transition temperature well 
in the strong-coupling regime. The IPT is a less accurate approximation for both the self 
energy and vertices for retarded interactions, but is much more accurate in determining 
transition temperatures because the errors in the self energy and the vertices tend to cancel 
in the determination of Tc. The the explanation for this behavior is that the IPT is an accu- 
rate approximation when third-order terms can be neglected. This occurs for the self energy 
of the Hubbard model at half filling, but is not true for the self energy when the interaction 
is retarded or the filling is different from 1, or for the irreducible vertex functions. Qualita- 
tively the IPT approximation is better than the same order conserving approximationEl for 
the self energy of the Hubbard model, but the transition temperatures are better approxi- 
mated (quantitatively) by the conserving approximation. For the Holstein model, the IPT 
is superior to the conserving approximation in all aspects except for the anomalous behavior 
in Tc near half filling and for large enough coupling strength. 
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FIGURES 

FIG. 1. Dyson equation for the self energy of the Holstein model. The solid lines denote the 

bare (electronic) Green's function and the wavy lines denote the phonon propagator. The self 
energy includes the Hartree and Fock contributions, the second-order dressing of the phonon line, 
the lowest-order vertex correction, and the second-order rainbow diagram (which arises from the 
self-energy insertion of the Fock diagram into the Fock diagram) . 

FIG. 2. The irreducible vertex functions in the CDW and SC channels. The CDW irreducible 
vertex function is shown in (a). Note that the vertex corrections (exchange diagrams) modify the 
interaction to lowest order in the CDW channel. The SC irreducible vertex function is shown in (b). 
The vertex corrections first enter at second order in the SC channel. In the IPT approximation the 
electron propagator is the bare Green's function Go, while in the IPT* approximation the electron 
propagator is the dressed Green's function G. 

FIG. 3. Comparison of the IPT (solid line) and IPT* (dashed line) to the QMC solutions (solid 
dots) for the Holstein model at half filling with phonon frequency = 0.5t*, interaction strength 
g = OAt*, and temperature T = t*/16. This example is generic for the weak-coupling limit. In (a) 
the self energy renormalization function Z{iujn) — 1 is plotted against the Matsubara frequency. In 
(b) the symmetric combination of the first column of the irreducible vertex function in the CDW 
channel is shown. Note that even at this weak a value of coupling strength the third-order diagrams 
must play an important role. 

FIG. 4. Comparison of the IPT (solid line) and IPT* (dashed line) to the QMC solutions 
(solid dots) for the Holstein model at half filling with phonon frequency J7 = 0.5t*, interaction 
strength g = Q.5t* , and temperature T = t* /9. This example is generic for the transition region to 
the strong-coupling limit. The self-energy renormalization function (a) and the irreducible vertex 
function in the CDW channel (b) are pictured. Note that in the limit where the strong-coupling 
effects begin to manifest themselves, the QMC vertex becomes increasingly attractive at low fre- 
quency transfer. 
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FIG. 5. Comparison of the IPT (solid line) and IPT* (dashed line) to the QMC solutions (solid 
dots) for the Holstein model at half filling with phonon frequency O = OM*, interaction strength 
g = 0.625t*, and temperature T = t* /Q. This example is generic for the strong-coupling limit. The 
self-energy renormalization function (a) and the irreducible vertex function in the CDW channel 
(b) are pictured. Note that in the strong-coupling regime the QMC vertex becomes quite attractive 
at low frequency transfer. 

FIG. 6. Transition temperature to the CDW-ordered state at half filling in the Holstein model 
at an intermediate phonon frequency {^l = OM*). The IPT approximation (solid line) is compared 
to the IPT* approximation (dashed line) and the QMC results (solid dots). Both approximations 
are similar to each other and are accurate at predicting the peak height and position. This occurs 
because the underestimation of the self energy is cancelling the underestimation of the vertex in 
the calculation of Tc. 

FIG. 7. Phase diagram of the Holstein model with = 0.5t* at three different coupling 
strengths {g = 0.4,0.5,0.625). The solid dots are the QMC solutions with CDW order, and the 
open triangles are the QMC results with SC order (the dotted lines are a guide to the eye). The 
kinks in the solid (IPT) and dashed (IPT*) lines occur at the CDW-SC phase boundaries. Note 
the anomalous behavior of the CDW transition temperature not being maximal at half filling sets 
into the approximate theories as the coupling strength increases. 

FIG. 8. Self energy diagrams for the Hubbard model thru fourth order. The solid lines denote 
the bare electronic Green's functions, and the dotted lines are the Coulomb interaction. The first 
two lines contain all of the bubbles and ladders thru fourth order. The middle two lines are 
the remaining diagrams that enter first at fourth-order. The last line includes the second-order 
self-energy insertions into the second-order diagrams. At half filling the odd-order contributions 
to S vanish, and each of the three fourth-order contributions on the same line yield the same 
contribution to S. 
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FIG. 9. Near cancellation of the fourth-order contributions to the self energy. The second-order 
contribution to the self energy (dashed line) is compared to each of the four fourth order contri- 
butions (dotted lines) and the net fourth-order contribution (solid line) at U = —t*, T = t*/20, 
and half filling. The dotted line denoted FLEX corresponds to the first line of the fourth-order 
contributions in Figure 8. The labels line 2 and line 3 denote the next two lines of fourth-order 
contributions in Figure 8, and the label IPT refers to the last line in Figure 8 corresponding to 
the second-order self-energy insertions into the second-order diagrams. Note that the total contri- 
bution of the fourth order diagrams is an order of magnitude smaller than any of the individual 
classes of diagrams in Figure 8, and that it is the same sign as the FLEX contributions. 

FIG. 10. Comparison of the IPT and IPT* approximations to the QMC solutions for the 
Hubbard model at half filling in the limit of weak-coupling {U = —t*, T = t*/20). The second-order 
(dashed line), third-order (dotted line), and fourth-order (solid line) approximations are compared 
the QMC solutions (solid dots). In (a) the self-energy renormalization function is plotted against 
Matsubara frequency. In (b) the even component of the first column of the irreducible vertex 
function in the CDW channel is plotted in the IPT approximation, and in (c) it is plotted in the 
IPT* approximation. 

FIG. 11. Same as in Figure 10, but with a stronger value of the couphng (U = —2t*, T = t*/8). 
The self energy is still reproduced accurately by all orders, but only the third-order IPT approxi- 
mation has the correct qualitative behavior for the vertex. 

FIG. 12. Same as in Figure 10, but with a stronger value of the coupling {U = —3t*, T = t*/6). 
The self energy is still reproduced accurately by all orders, but only the third-order IPT approxi- 
mation has the correct qualitative behavior for the vertex. 
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FIG. 13. Transition temperature to the CDW-ordered (and SC-ordered) state in the Hubbard 
model at half filling. The second-order (solid line), third-order (dotted line), and fourth-order 
(solid line) IPT (a) and IPT* (b) approximations are compared to the QMC results (solid dots). 
Note that the transition temperature curves all have a turnover (except for the third-order IPT 
approximation) because the vertex is underestimated as the coupling strength increases. The 
third-order IPT approximation does not turnover because the vertex is overestimated. 

FIG. 14. Transition temperature to the SC state in the Hubbard model as a function of electron 
concentration at U = —1. The second-order IPT (solid line), third-order IPT (dotted line), sec- 
ond-order IPT* (dashed line), and third-order IPT* (chain-dotted line) are compared to the QMG 
results (open triangle). Note that both second-order approximations display anomalous behavior 
near half filling even for this small a value of the interaction strength. 
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